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the couplings of the effective theory and the parameters N T ,f3 of the original thermal lat- 
tice theory, which can be systematically improved. We then investigate the deconfinement 
transition for the cases SU{2) and SU (3) by means of Monte Carlo simulations of the effec- 
tive theory. Our effective models correctly reproduce second order 3d Ising and first order 
phase transitions, respectively. Furthermore, we calculate the critical couplings (3 C {N T ) and 
find agreement with results from simulations of the 4d theory at the few percent level for 
N T = 4 - 16. 
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1 Introduction 

Non-abelian gauge theories at finite temperature are inherently non-perturbative, due to 
infrared problems for soft gauge fields [1]. The dynamical appearance of different scales, 
ttT, gT, g 2 T with gauge coupling g, has motivated effective theory methods, most notably 
dimensional reduction [2, 3]. The idea is to integrate over the hard scales perturbatively, 
whereupon a three-dimensional effective theory arises which is then easier to solve by 
non-perturbative means. A particularly successful application is the investigation of the 
electroweak phase transition as a function of the Higgs mass using the 3d effective theories 
[4, 5], which produced results at the few percent accuracy level compared to simulations of 
the full 4d theory [6]. 

In QCD, the same technique can be applied to observables in the deconfined phase only, 
because the reduction loses validity as the deconfinement transition is approached. An in- 
vestigation of the phase diagram by similar methods would be particularly desirable, since 
lattice Monte Carlo studies at finite baryon density are beset by the sign-problem. Unfor- 
tunately, the QCD phase structure is not inherited by the dimensionally reduced model, 
which loses the Z(N) symmetry of Yang-Mills theory in the perturbative reduction step 
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[7]. This has motivated alternative approaches, where the most general Z(A r )-symmetric 
theory is written down and the couplings are then matched to those of the full theory 
by calculations of particular observables [8, 9]. While a succesful description of the phase 
transition for SU(2) can be achieved in this way [10] (see also [11]), in the case of SU(3) 
there remain so far several open couplings that cannot be matched easily. 

In this paper, we study the possibility of reducing the full theory to an effective theory 
by using strong coupling expansions on a lattice. Such an approach was conjectured to be 
sensible in [12] and has been explored earlier in the literature [13-18]. It results in a three- 
dimensional effective theory of Polyakov loops. A common simplification was the neglect 
of spatial plaquettes, which has been argued not to influence the universal behaviour of 
the theory. Work where authors went beyond the spatial strong coupling limit is [19], see 
also [20] for recent developments on Polyakov loop extended strong coupling lattice QCD 
with staggered fermions. 

Here we significantly extend this approach and calculate longer strong coupling series 
for the effective couplings, which are thus valid beyond the spatial strong coupling limit. 
Strong coupling series have a finite radius of convergence, which in this case is given by 
the deconfinement transition. Hence, our effective theory is valid in the confined phase and 
complementary to weak coupling approaches. We also investigate the influence of next-to- 
nearest neighbour Polyakov loop interaction terms which arise in higher orders. For our 
first studies we consider the pure gauge theories SU(2) and SU(3). After derivation of the 
effective theories of complex scalar fields in 3d, corresponding to the traces of the Polyakov 
loop and its hermitian conjugate, we study them by means of Monte Carlo simulations. 
In particular we determine the order and the critical couplings of the deconfinement phase 
transition and relate them to those of the original theories. We find qualitative and quan- 
titative agreement with less than 6% deviations for the (3 C (N T ) up to N T = 16, suggesting 
that these models are useful for the study of continuum physics. 

The paper is structured as follows: in section 2 the effective actions for the two gauge 
groups are derived, with detailed discussion about the higher-order contributions to the 
effective couplings; then, in section 3, we present the Monte Carlo implementation of 
the effective models together with numerical results; section 4 provides the corresponding 
predictions for the original 4d theories and in section 5 we conclude. 

2 Derivation of the effective theory 
2.1 General strategy and SU{2) 

Consider the partition function of a (3 + l)-dimensional lattice gauge field theory at finite 
temperature IT = t4^J with gauge group SU(N) and Wilson's gauge action 



Ip J 

Finite temperature and the bosonic nature of the degrees of freedom imply the use of 
periodic boundary conditions in the time direction. 
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In order to arrive at an effective three-dimensional theory, we integrate out the spatial 
degrees of freedom and get schematically [16]: 



Z 



[dU ] exp [-S'eflf] 



-S cS = In J [dUi] exp A Y, (tr U p + tr E/j) 



= X\Si + X2S2 + . . . . 



(2.2) 



We expand around /3 = and arrange the effective couplings A n = X n ((3,N T ) in increasing 
order in f3 of their leading terms. Thus, the X n become less important the higher n. As we 
shall see, the interaction terms S n depend only on Polyakov loops 



L 4 



tiWj 



tr H U ( 



x 



r) 



(2.3) 



This is the reason for a "dimensional reduction" occurring here, as the time dimension is 
now implicit in the variables of the effective theory, which are fields defined on the spatial 
lattice. With sufficiently accurate knowledge of the relations \ n (/3,N T ), we are able to 
convert the couplings of the three-dimensional theory to those of the full theory. In this 
work we are mainly interested in the deconfinement transition. Determining the critical 
parameters A njC of the effective theory then gives a whole array of critical f3 c (N T ) for - in 
principle - all N T . In the following we calculate strong coupling, i.e. small /3, expansions of 
the leading X n . 

Since the calculations are quite similar for different numbers of colours, we now spe- 
cialise our derivation to the simpler case of SU(2) and later provide the necessary changes 
for 577(3). Using the character expansion as described e.g. in [21, 22], the effective action 
according to Eq. (2.2) can be written as 



S e ff = In 



/ Wi] n 



1 + )d r a r (P)xr(U P 



(2.4) 



where the sum extends over all irreducible representations r with dimension d r and char- 
acter Xr- The expansion coefficients o r (/3) are accurately known [21] and in the following 
we use u = chf as expansion parameter instead of f3 for its better apparent convergence. 
The logarithm in this definition allows us to use the method of moments and cumulants 
[23], and we get the following cluster expansion 



S c ff 



E ■ 

c=(xl H ) 



Hi 



(2.5) 



[dUi] J| d rp a rp Xr p (U p ) , 



where the combinatorial factor a(C) equals 1 for a single polymer Xj and —1 for two non- 
identical connected polymers. For clusters consisting of more than two polymers, a(C) 
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Figure 1. First graph with a nontrivial contribution after spatial integration for a lattice with 
temporal extent N T = 4. Four plaquettes in the fundamental representation lead to an interaction 
term involving two adjacent fundamental Polyakov loops Li and Lj. 

depends on how these polymers are connected. In contrast to earlier calculations [24, 25], 
where also the temporal links were integrated over, the contribution $ still depends on 
Wilson line variables Wj, which results in the interaction terms for the effective action, 
Eq. (2.2). Our task is then to group together all graphs yielding the same interaction 
terms up to some order in /3, and this finally gives the strong coupling expansion of the 
corresponding effective coupling A n . 

2.2 The spatial strong coupling limit 

In earlier strong coupling calculations at finite temperature it has been customary to neglect 
spatial plaquettes. Let us briefly discuss this limit as it gives some important insights for 
the following. Neglecting spatial plaquettes, the integrations in Eq. (2.4) can be done 
applying the group integration rules 



where we replaced the integration over all temporal links with an integration over the 
Wilson line variables Wi as in [26]. This is justified since the integrand depends on temporal 
link variables only through the Wj. 

Some important observations in this limit are: 

• The summation extends only over pairs of nearest neighbours <ij>, i.e. next-to- 
nearest neighbour interactions vanish without the inclusion of spatial plaquettes. 

• Contributions from higher representations start with higher powers of f3. 

• The exponential function has cancelled, hence if we insist on having a partition 
function of the form Z = j exp[— S], we have to introduce a logarithm in the action. 

2.3 Leading order effective action 

The leading order result of the effective action has first been calculated in [13] and corre- 
sponds to a sequence of iV T plaquettes that wind around the lattice in temporal direction, 
cf. figure 1. Its contribution is given by: 






on each spatial link and we get the partition function [14]: 







<ij> 
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Figure 2. Left: First correction to the leading order graph, proportional to N T . Middle, right: 
Repetitions of this decoration. 



Hence, to leading order the first coupling of the effective theory is \\(u, N T ) = u Nt . 

For additional terms of the series for Ai, we can use most of the graphs that also appear 
in the strong coupling expansion of the Polyakov loop susceptibility [25] . These corrections 
involve additional plaquettes, are hence of higher order in u and we call these attached 
plaquettes decorations. Let us note that repetitions of lower order decorations attached to 
planar graphs exponentiate and we can write 



Ai(it, N T ) = u Nt exp \n t P(u, N t ) 



(2.10) 



with some polynomial P(u,N T ). This can be seen e.g. from the graphs shown in figure 2 
and their corresponding contributions: 



Left: 
Middle: 

Right: 



$3 



II 



N T 



11 



N T 



4lN t u St ; 

| (47V T n 4 ) • 4(iV r - 3)i 



Si 



,Nr 



4iV r ii 4 • 3N r u 4 



Si . 



Combining the three parts, we can write this, up to higher orders, as 



$1 + $2 + $5 



u exp 



N T (4n 4 - I2u 8 ) 



Si 



(2.11) 



(2.12) 



Of course, the polynomial in the exponential is only part of the complete result for Ai to 
that order. For example, there are other graphs contributing to order u e which are still 
missing in this correction. Nevertheless, in this exponentiated form the effective coupling 
corresponds to a partial resummation of higher order terms which may be expected to 
improve convergence behaviour. Let us remark that such an exponentiation has been 
observed also for the strong coupling expansion of the string tension. 

Carrying out the calculations, we get the following results through order u 10 in the 
corrections relative to the leading order graph: 

49044 



Ai(it, 2) = u 2 exp 
Ai(w, 3) = -u 3 exp 
Ai(w, 4) = u 4 exp 
Xi(u,N T > 5) 



N T 

u T exp 
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N T ( 4u A 



u 134 „ 
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(2.13) 



For smaller iV T some graphs do not contribute since the temporal extent of their decoration 
is > N T so that they do not fit into the lattice. The coefficients of the order n in the 
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Figure 3. Shape of the most elongated graph contributing to Ai at order n = 2(q + 1), here in the 
case q = 3: only systems with N T > q can actually accommodate for it. 

exponents of the effective couplings, Eq. (2.13), reach their N T = oo values as soon as 
N T = n/2. This can be understood as follows. Among all graphs contributing to the 
coefficient at order u n , the most elongated in the time direction is of the type depicted in 
figure 3, with q consecutive lifted plaquettes, which brings inn = 2(q+l) additional powers 
of u. Such a graph is included only if N T > q + 1 (in the moment-cumulant formalism a 
given polymer cannot occupy twice the same plaquette). 

2.4 Higher order terms 

There occur several types of higher order graphs: larger numbers of loops involved, Polyakov 
loops at distances larger than one and Polyakov loops in higher dimensional representations. 

We begin by considering powers of the leading order term. Inspection of higher order 
terms shows that one can arrange a subclass of these terms in the following manner 

£ (\\LiLj - |LfL| + Ml^s = £ (i + X lLlLj ) . (2.14) 

<ij> ^ ' <ij> 

Thus, there are graphs that reproduce the emergence of the logarithm just as in the spatial 
strong coupling result, Eq. (2.8). In contrast to that case we have now the full effective 
coupling Xi(u,N T ) appearing in the logarithm instead of only its leading order term u T , 
which results if we restrict Eq. (2.8) to the fundamental representation. To see this, one 
calculates the corresponding graphs with L 2 L 2 or LfL^, and the combinatorial factor a(C) 
of Eq. (2.6) gives the correct prefactors for the series to represent a logarithm. 

Next, let us consider couplings pertaining to next-to-nearest neighbour interactions. 
These appear once additional plaquettes are taken into account. Naively, the leading 
contribution should correspond to a planar graph with Polyakov loops at distance two. 
However, this graph is precisely cancelled by the contribution of the nearest-neighbour 
graph squared and its associated combinatorial factor —1 (figure 4). The leading non-zero 
contribution therefore comes from L-shaped graphs and is given by 

A 2 («, N T )S 2 = N T (N T - l)u 2N ^ +2 L k Li , (2.15) 

[kl] 

where we have two additional spatial plaquettes (figure 5) and we sum over all pairs of 
loops with a diagonal distance of y/2a, abbreviated by [kl]. With the same steps leading 
to Eq. (2.14), we finally arrive at the SU(2) partition function 



Z= [dW] J] [1 + \\LiLj] H [1 + \ 2 L k U] 

<ij> [kl] 



(2.16) 



Figure 4. Two graphs without spatial plaquettes that cancel against each other due to different 
signs. 




Figure 5. Leading order graph contributing to the interaction between loops at a distance \/2a, 
denoted with [kl] in the summations. 



Let us further point out that corrections to the planar next-to-nearest neighbours, having 
a distance of 2a, give to leading order 

X 3 (u,N T )S 3 ^u 2N ^ +6 Yl L - L "' ( 2 - 17 ) 

«mn>> 

with an obvious notation for straight line next-to-nearest neighbours. The leading order 
graph is given in figure 6. Since these corrections are of higher order in u than the ones 
with distance \/2a, we will omit them in this work. As we shall see later, the effective 
theory in case of SU{2) works quite well even without next-to-nearest neighbours, hence 
we will quote explicit results for the coupling A2 for the case of SU(3) only. 

Finally, we include some remarks about the Polyakov loops in higher dimensional 
representations. Consider, e.g., the adjoint Polyakov loop: the leading order term emerging 
from a strong coupling expansion is 

2 2 16 

XaS a = £ XamXaiWj) , V = -U 2 + -U 4 + — U 6 + . . . 

<ij> 

and hence A a ~ u 2Nt , which is formally of lower order than the coupling A2 . To next-to- 
leading order (valid for all N T > 2) we have 

X a = v N ^(l + N T l- + ..). (2.18) 



3 v 

Effects of higher representations have also been investigated in the literature [17, 18, 27]. 
2.5 The effective action for SU(3) 

In the case of SU (3) the same steps as for SU(2) apply. The only difference we have to 
keep in mind is that SU(3) also has an anti-fundamental representation and consequently 
there is also a complex conjugate Polyakov loop variable L*. Thus we get the one-coupling 
and two-coupling partition functions 

Zt = j [dW] H [l + Ai (LiL* + L*Lj)] , (2.19) 

<ij> 

Z 2 = j [dW] n [1 + Ai {UL* + L*Lj)] Yl [1 + A 2 {L k L\ + V k U)] . (2.20) 

<ij> [kl] 
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Figure 6. Leading order graph of next-to-nearest neighbours with distance 2a. 



The effective coupling Xi(u, N T ) is obtained as (for this gauge group, in view of the 
comparison that will be made with available 4d Yang-Mills simulation data, we consider 
only even values of N T ; however, analogous formulae hold for odd N T ): 

Ai(2, it) = u 2 exp [2 (4u 4 + 12u 5 - 18u 6 - 36u 7 
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2 10 
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Ai(JV T > 6, it) 
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295 
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1851 9 1055797 



-u a + 



5120 
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(2.21) 



For the first terms of the next-to-nearest neighbour coupling \2(N T ,u) we find 

A 2 (2,n) = u A 



2u 2 + 6« 4 + 31u e 



u 



12u 2 + 26n 4 + 364u 6 



.12 



U 



30n 2 + 66u 4 



A 2 (4,u 
A 2 (6, u 

X 2 (N T > 8,u) = u 2Nt [N t (N t - l)u 2 ] , 
while the leading coupling of adjoint loops is (valid for iV T > 2) 



2 v 



81 



-u H u + 

S 32 



(2.22) 



(2.23) 



3 Numerical simulation of the effective theories 



In this section we proceed to evaluate the effective theories by Monte Carlo methods. 
After introducing a suitable algorithm, we discuss the strategy employed to extract critical 
couplings and the properties of phase transitions, and present detailed numerical results. 

3.1 The one coupling model 

For the purpose of numerical simulations, a further simplification is achieved by using the 
trace of the Polyakov loops for the path integral measure as degrees of freedom (complex 
numbers instead of matrices), and rewrite the one-coupling partition function for SU(3), 
Eq. (2.19), 

Z = ( II / dL *) e ~ ScB ; S ^s = - Yl Ml + 2AiReLjL*) - . (3.1) 

x <ij> x 
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Figure 7. Domain for the complex numbers L x in the SU (3) effective theory. 



The complex numbers L x represent traces of SU(3) matrices, and thus \L X \ < 3. The 
complex domain of L x is thus the "three-pointed star" with radius 3, figure 7. 

The potential term in S e g is the Jacobian induced by the Haar measure of the origi- 
nal group integration and its actual form depends on the chosen parametrisation for L x . 
Following [28], we write the trace of a SU(3) matrix by rotating it to its diagonal form, 

L x (0, 4>) = e ie + e** + e~ i( - e+ ® , -it < 9, <p < +ir . (3.2) 
In this case the potential is 

V x = i log(27 - 18IL,! 2 + SRbLI - \L X \ A ) . (3.3) 
The integration measure actually used in our simulation then takes the form 

[ dL x e v * = f + \(k x I ' " &0 x e v * . (3.4) 

J J— 7T J— 7V 

When working on the 577(2) theory, —2 < L x < +2 is a real number and we simply have 

dL x e v " ,V* = \ log(4 - Ll) . (3.5) 
3.2 A "sign problem" and its solution 

The effective theory described so far has lower dimensionality and simpler degrees of free- 
dom per site than the initial Wilson action, suggesting a straightforward local Metropolis 
update based on accept/reject steps. This ideal situation is somewhat complicated by some 
sort of a "sign problem" . 

Consider the effective partition function in Eq. (2.19). For couplings larger than some 
threshold, A > A T , there are gauge configurations yielding negative contributions to Z , 
or, equivalently, for which the logarithms appearing in the effective action have negative 
arguments. Although Z remains overall positive, the Metropolis update ceases working 
when approaching this threshold value of the coupling. 

For the SU(3) action ReLjL* can go down to —9/2 and the threshold coupling is 
X T = 1/9, while for SU(2) we have X T = 1/4. In this work our interest is in the phase 
transition, and we shall see that for SU(3) the critical coupling Ai )C happens to be close to 
X T . This sign problem is then a practical concern and a solution is called for. Fortunately, 
in the case of SU(2) we have Ai jC ~ 0.2 < A T = 1/4 and this problem can be ignored. We 
exploit this fact to test our workaround for the problem against the full solution. 
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Figure 8. Left: Distribution of L for small and large Ai on a lattice with N s = 6 and M = 1. 
Middle, Right: Expectation value of \L\ and its susceptibility. The vertical line marks the infinite- 
volume transition. 



Our approach to overcome this problem is the following: we Taylor-expand the loga- 
rithm in the effective action to some order M in powers of q = AiReLjL*, which effectively 
undoes the resummation as in Eq. (2.14), 

*2° = - E v. - E - 2 « 2 + ^ 3 - ^ + f - • • • - ^ M w qM ) ■ (3 - 6) 

x <ij> 

We thus obtain a family of models which are, by construction, free of the mentioned 
problem, and are expected to converge to the full theory as M — > oo. The strategy is then 
to find the series of critical points Ai jC (M) and to look for their convergence behaviour as 
M increases. 

In order to have a critical model, moreover, one must truncate the expansion at an 
odd power M: for even truncations, indeed, the slope of the effective action as a function 
of q, corresponding to a force, vanishes for 2q < 1, and never triggers a symmetry breaking 
by favouring aligned neighbours. Let us observe, incidentally, that this limitation is absent 
in the SU (2) case, exactly because the transition region is far from the threshold coupling 
and the value of \\LiLj is always too far from 1 for this phenomenon to take place. 

3.3 Phase structure, critical coupling and finite size analysis 

Our first task is to establish the phase structure of the effective theory, where we focus 
on the physically interesting case of SU(3). Based on the global Z{2>) symmetry of the 
model, one expects spontaneous breaking of that symmetry for some critical value of the 
coupling Ai jC . Figure 8 shows the behaviour of the field variable L as a function of Ai. As 
expected from the 4d parent theory, there is indeed a transition from a disordered or mixed 
phase, with values of L scattering about zero, to an ordered phase at large coupling where 
the three Z(3)-phases are populated separately. In the thermodynamic limit, one of these 
vacua will be chosen and the symmetry is broken spontaneously, (L) = for Ai < Ai )C 
and (L) ^ for Ai > Ai jC . Correspondingly, the expectation value of \L\ rises abruptly 
at some critical coupling Ai jC , as shown in figure 8 (middle). On a finite size lattice, the 
phase transition is smoothed out, non-analyticities are approached gradually with growing 
volume, as the figure illustrates. 
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The general technique to locate the infinite- volume critical coupling, Ai )C , is based on 
a finite-size scaling analysis. A variety of cubic systems are simulated and for each one a 
pseudo-critical point Ai jC (iV s ) is found. As a function of lattice side length N s , close enough 
to the thermodynamic limit these scale as 

Ai, c (JV.) = A liC + bN-V» . (3.7) 

The relevant universal values for the critical index are v = 0.63002 for 3d Ising [29] and 
v = 1/3 for a first order transition. In practice we found it sufficient to use spatial volumes 
N s = 6, . . . , 12 in order to extrapolate to the thermodynamic limit. In this way, the whole 
data production can be carried out in a handful of days on a modest desktop PC. 

The pseudo-critical coupling can be defined in a variety of ways: one possibility is to 
employ the energy observable: 

E = -^S' eS , (3.8) 

where the prime denotes omission of the potential term V x . For the M = 1 model, this 
coincides with the usual energy, 1 while at higher truncations (as well as in the non-truncated 
formulation) S e g is non-linear in the coupling. For this reason, the energy is used only in 
the M = 1 particular case. Other natural observables are based on the modulus of the 
Polyakov loop, \L\. We consider the Binder cumulant 

« = 1-^, (3-9) 

whose minimum defines a pseudo-critical coupling, and the susceptibility 

X (\L\) = ({\L\ - (\L\)Y) , (3.10) 

whose maximum is taken as pseudo-critical coupling, cf. figure 8 (right). Having various 
definitions of pseudo-critical couplings, one can check for the mutual consistency of the 
infinite-volume critical points coming from them. 

3.4 Critical coupling and order of the transition for SU (3) 

The truncated theories with M = 1,3,5 were simulated on lattices with spatial sizes 
N s = 6,8,10 (plus N s = 12 for the M = 1 theory). For each volume, ~ 30 values of 
the couplings are sampled by ~ 10 6 update sweeps each. Measurements were taken every 
~ 30 updates. 

We begin our presentation of results with the SU(3) family of effective theories, where 
we find a first-order transition regardless of the particular truncation employed. A practical 
difficulty is the related occurrence of metastabilities with extremely long thermalisation 
times oc exp(a/V;? ). The local Metropolis update requires a number of iterations akin to a 
tunnelling time in order to dispel configurations with two halves of the system associated 
with different vacua. To take care of this problem, very long trajectories are inspected 
by cutting out an initial block of configurations until stability is found in the resulting 
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Figure 9. Left: Position of the minimum of the Binder cumulant B(E) for SU(3),M = 1, for 
different lattice sizes. The horizontal line is the thermodynamic limit resulting from the fit to 
Eq. (3.7). Right: Behaviour of B m i n (N s ), along with its thermodynamic limit obtained with the 
0(N^r 3 ) scaling law and the independent estimate B^ from the \L\ histogram. Also the second- 
order limit value 2/3 is shown. 



averages: for instance, a system with size N s = 16 would require, around criticality, ~ 10 
update sweeps to thermalise. 

First we consider the model with M = 1. The first-order nature of the transition 
is established by fitting the pseudo-critical couplings to the scaling law, Eq. (3.7), with 
v = 1/3, as shown in figure 9 (left). The same conclusion is reached when analysing the 
scaling of the minimum i? m ; n (iV s ) of B(\L\), figure 9 (right): this quantity, in the N s — > oo 
limit, approaches the saturation value 2/3 if and only if the transition is second-order, 
while in case of a first-order point the thermodynamic limit is given by [30]: 

Boo = ---(W-r¥)^ (3-ii) 

3 12V|L| 2 \L\xJ ' v ' 

where |L|i and \L\i are the local maxima of the \L\ double-peaked histogram. Inspection of 
the \L\ distribution shows indeed a two-peaked shape (figure 10) that becomes narrower and 
higher at increasing system volumes. It is then possible to give an independent estimate 
for Boo and compare it with that from the scaling of the minimum. The comparison 

(figure 9, right) confirms indeed the first-order nature of the transition, while some slight 
mismatch (less than two standard deviations) remains between the two estimates. This is 
most likely due to the presence of subleading terms in the scaling of B m [ n . As demonstrated 
in [31], the scaling law to NLO has the form 

B min (N s ) = £00 + B^N- 3 + B^N~ 6 , (3.12) 

and it is necessary to explore very large lattices to observe the leading-order scaling. 

1 Also in the A — > limit, at any truncation, this definition recovers the standard energy term. 
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Figure 10. Left: Histogram for \L\ of the SU(3) theory on N s = 10 with M = 1. The three curves 
have been measured with Ai slightly below, exactly at and slightly above pseudo-criticality. The 
horizontal lines are guides to the eye and mark the peak on the right as being, for the pseudo-critical 
curve, exactly three times as high as the one on the left. Right: Histogram of the SU(3) \L\ for 
four different system sides N s at the true critical point. The size-dependence of the peaks' position 
was taken into account in estimating with Eq. (3.11). 
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Figure 11. Critical points in the thermodynamical limit (horizontal lines) and pseudocritical values 
(data points) on three spatial sizes, for the SU{3) one-coupling model truncated at M — 1, 3, 5. 

In the next step we need to investigate the behaviour of the models with higher M. 
Again we observe first order transitions, which become sharper with increasing M. More- 
over, finite-size effects are stronger for higher M, figure 11. The critical couplings identified 
for the M = 1, 3, 5 effective theories in the thermodynamic limit are also quoted in figure 
11. Judging from these three values, the series seems to be rapidly converging, with only 
~ 3% difference between M = 3,5. The residual difference between this estimate and the 
M = oo critical coupling is completely subdominant compared to the other systematic 
errors contributing to the final results. Also, the direct comparison with the SU(2) case 
below, where the M = oo data are directly available, supports a rapid convergence, figure 
13. 
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Figure 12. SU(2) with M = 1; Left: Location of maximum of the -E-susceptibility, fitted to the 
3D Ising scaling law. Right: Histogram for \L\ at N s = 20. 



3.5 Critical coupling and order of the transition for SU(2) 

In the case of the SU{2) models, the same analysis was carried out. Since here we find 
second-order transitions, the metastability problems are absent and the thermalisation 
times are orders of magnitude shorter (e.g. around 4000 update sweeps at N s = 16 around 
criticality). This enabled us to study larger system sizes up to N s = 28 without much 
computational effort. 

The same approaches as before clearly exhibit the second order nature of the transition, 
using the scaling law for the pseudo-critical points and the histogram for \L\ (see figure 
12), as well as the corresponding Binder cumulant thermodynamic limit which in this case 
reaches smoothly the value 2/3. Figure 12 shows the results for the M = 1 model. However, 
as indicated before, in the case of SU (2) we can directly work with the full M = oo model 
and thus check for the systematic errors when a finite M truncation is used. This is shown 
in figure 13, explicitly we find 

Ai, c (M = 1) = 0.195374(42) ; A 1]C (M = oo) = 0.21423(70). (3.13) 

Indeed we observe rapid convergence on the M = oo results, in accord with our findings for 
SU(3). We thus conclude that the truncation in M does not affect the results significantly 
and proves to be a viable way of dealing with the numerical difficulties discussed above. 

3.6 Two-coupling models for SU (3) 

In this section we study the influence of including a second coupling. We start with 
the leading next-to-nearest-neighbour interaction, specified by the model with a second 
coupling A2, Eq. (2.16). For SU(S), the two-coupling partition function reads 

Z = ( II / dL -) II 0- + 2AiReL i L*) J](l + 2X 2 ReL k Lf)e^ v * . (3.14) 

* <ij> [kl] 

As discussed previously, we neglect a third coupling related to next-to-nearest-neighbour 
interaction between loops at distance 2a since it starts at higher order than A2. 



-14- 



data from B{ L ) i — 

data from xi L ) L ~ 
M=infinity, B L) 
M=infinity, %{ L ) 

, : ■ ; • 



1 2 3 4 5 6 7 8 

M (truncation order) 

Figure 13. Xi^ c (N s = 8) for SU(2) and various truncations. Estimates from j5(|L|), x(|L|) are 
shown; horizontal lines mark the values for the model without truncations (M = oo). 

In a natural extension of the critical point of the one-coupling model, there is now a crit- 
ical line Ai )C (A2) separating the ordered from the disordered phase in the two-dimensional 
parameter space. However, not all points on this line are related to the physics of the 4d 
thermal theory. Once a particular N T is fixed, both Aj depend on the expansion parameter 
u{(3) alone. Eliminating u, the curves can be constructed, which represent the 

parameter space describing the physics of the 4d thermal model for a given N T . The critical 
couplings relevant for us are thus the intersections between the critical line of the effective 
two-coupling model and the curves specifying the map to a particular N T lattice. 

We remark that the second factor in the partition function, Eq. (3.14), has the same 
kind of "sign problem" as the first one, calling for an analogous solution. Some attention 
needs to be paid to truncating the expansion of the two logarithms in a consistent way. 
We include all terms which contribute to the desired power of u from both expansions, for 
all values of N T . Based on our previous results, we want to keep terms up to M = 3 for 
the Ai part, and thus consider the leading second coupling term M = 1, which we denote 
with {M\,M2) = (3,1). Corrections by higher powers M2 are formally of higher order in 
u, and we have checked numerically that their effect is completely negligible. 

In order to locate the critical line Ai )C (A2), we fix 11 different values of A2 in the range 
[0.002, 0.007] and scan in Ai in order to identify the pseudo-critical points, which can then 
be extrapolated to the thermodynamic limit. It turns out that, on a given volume, Ai iC (A2) 
is well described by a linear interpolation, cf. figure 14 (left). Moreover, the inclusion of 
the second coupling appears to diminish finite size effects. Since all three volumes shown 
in figure 14 (left) are consistent within errors, we quote the result from the largest lattice 
for the critical line, 

A ljC = a + 6A 2 with a= 0.10628(8), b = -1.891(4) . (3.15) 

The value of a is consistent with the M = 3 critical point in the one-coupling theory. 

Also shown in figure 14 (left) are the lines \\ Nt \\2) for N T = 4, 6, 8. One observes that 
they rapidly accumulate towards A2 — 0. Only for the lowest values of N T , corresponding 
to coarse lattices in the 4d theory, does the adoption of the two-coupling model make 
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Figure 14. Critical line in the SU(3) two-coupling space, determined from x(|L|). Dashed lines 
give the parameter space representing a 4d theory with fixed N T . Left: (Ai, A2). Right: (Ai, A a ). 



any difference. For finer lattices, of interest for continuum physics, the results are within 
statistical errors indistinguishable from the the simpler one-coupling theory. 

Next, we are considering a two coupling theory with Ai, A a . In this case the partition 
function reads 

Z = ( II / dL *) II ( X + 2AiReL,L*) [] e A a (TrW^)(TrW^) e E, v x (3 16) 

x <ij> <ij> 

Here Tr*- a ^ denotes the trace in the adjoint representation, Tr^W = \TrW\ 2 — 1, which can 
be used to rewrite the action in terms of fundamental loops. For the numerical evaluation, 
we again expand the Ai,A a terms to (Mi = 3, Mi = 1), and proceed in complete analogy 
as in the case discussed above. The result for the critical line in this two-coupling space is 
shown in figure 14 (right). Here we find 

A liC = a + 6A a with a= 0.10637(15), b = -1.422(22) . (3.17) 

Once more, the set of curves intersecting the critical line correspond to lines of fixed iV T in 
the 4d theory. We observe that A a has slightly larger effect than A2 at fixed N T , in accord 
with the fact that it starts at lower order in u, cf. Eq. (2.23). Nevertheless, its influence 
is smaller than that of the strong coupling truncation in Ai, as we shall see, and hence 
negligible at this order. 

4 Mapping back to 4d Yang-Mills 

Having established the critical couplings for our effective theories and tested their reliability, 
we are now ready to map them back to the original thermal Yang-Mills theories by using 
Eqs. (2.13, 2.21). In Tables 1, 2 we collect the values for the critical gauge couplings, /3 C , 
obtained in this way from the effective theories and compare them to the values obtained 
from simulations of the full 4d theories for SU(2), SU(3), respectively. 

The agreement is remarkable in all cases, with the relative error of the effective theory 
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AT 


M = 1 


M = OO 


4d YM 


3 


2.15537(89) 


2.1929(13) 


2.1768(30) 


4 


2.28700(55) 


2.3102(08) 


2.2991(02; 


5 


2.36758(40) 


2.3847(06) 


2.3726(45) 


6 


2.41629(32) 


2.4297(05) 


2.4265(30) 


8 


2.47419(22) 


2.4836(03) 


2.5104(02) 


12 


2.52821(14) 


2.5341(02) 


2.6355(10) 


16 


2.55390(10) 


2.5582(02) 


2.7310(20) 



Table 1. Critical couplings (3 C for SU{2) from two effective theories compared to simulations of 
the 4d theory [32-34]). 



N T 


M = 1 


M = 3 


M 1 ,M 2 (X 2 ) = 3,1 


Mi,M 2 (A a ) = 3,l 


4d YM 


4 


5.768 


5.830 


5.813 


5.773 


5.6925(002) 


6 


6.139 


6.173 


6.172 


6.164 


5.8941(005) 


8 


6.300 


6.324 


6.324 


6.322 


6.0010(250) 


10 


6.390 


6.408 


6.408 


6.408 


6.1600(070) 


12 


6.448 


6.462 


6.462 


6.462 


6.2680(120) 


14 


6.488 


6.500 


6.500 


6.500 


6.3830(100) 


16 


6.517 


6.528 


6.528 


6.528 


6.4500(500) 



Table 2. Critical couplings /3 C for SU(3) from different effective theories compared to simulations 
of the 4d theory [33, 35]). 

results compared to the full ones shown in figure 15. The comparison of alternative trun- 
cations of the logarithm shows once more that it has almost no influence on the accuracy 
of the final result, as described earlier. Interestingly, there appears to be a 'region of best 
agreement', with the deviation growing both for small and large N T . We ascribe this to 
the fact that there are two competing systematic errors, as discussed earlier: the validity 
of the strong coupling series for a given coupling Aj is better the smaller (3 and hence N T , 
whereas the truncation of the next-to-nearest neighbour interactions gains validity with 
growing N T . In particular in the case of SU(3), there appears to be a cancellation of the 
two kinds of systematics, rendering the effective description better for the original theory 
on finer lattices. 

In order to asses the systematics of the strong coupling series, figure 16 (left) shows 
the difference in (3 C based on the one-coupling model when series of different depth are 
used for Ai(/3). Satisfactory convergence behaviour is observed. For the case of SU(2) a 
comparison with a non-perturbative derivation of the effective theory is possible. Figure 16 
(right) compares Ai )C in terms of our strong coupling series with a determination by inverse 
Monte Carlo [17] (see [18] for SU(3)). Again, we observe good convergence of the strong 
coupling series to the full result. However, these plots also illustrate that the error due to 
truncation of the strong coupling series of the first coupling is larger than the neglect of 
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Figure 15. Relative error of f3 c predicted by the effective theories when compared to simulations 
of the 4d theories, for SU{2) (left) and SU(3) (right). 
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Figure 16. Left: (i c for 5/7(3) from Ai, c (/3) evaluated to 8th, 9th and 10th order, Eqs. (2.21). 
Right: Ai(/3) for SU(2) and N T = 4 from strong-coupling expansion and inverse Monte Carlo [17]. 



higher couplings. We also attempted to improve the convergence of the series with Pade 
analysis, which however was not particularly successful, apart from confirming the order 
of magnitude of the estimated systematic uncertainties. 

Figure 16 (right) illustrates the range of validity of the strong coupling derivation of the 
effective couplings. The non-perturbatively determined Ai(/3) appears to change curvature 
at j3 c , whereas the estimates based on the strong coupling series do not. This is consistent 
with /3 C marking the radius of convergence also for the series expansion of the effective 
coupling Ai. Thus, the inverse Monte Carlo approach has a wider range of validity whereas 
the series approach furnishes analytically known mappings between the full and effective 
theories. 

5 Conclusions 

Employing strong coupling expansions, we have derived a dimensionally reduced, centre- 
symmetric effective description for lattice pure gauge theories at finite temperature. The 
effective theory is formulated in terms of scalar Polyakov loop variables and does not involve 
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any matrix degrees of freedom. Contrary to earlier derivations based on the lattice theory, 
we have included spatial plaquettes systematically. The effective theory has an infinite 
number of interaction terms with increasing degree of non-locality. Each of its couplings 
can be computed as a series in the original lattice gauge coupling (3 for any value of N T . We 
have explicitly calculated up to seven terms in the nearest-neighbour interaction and up 
to three terms in the next-to-nearest neighbour interactions. Furthermore, we have shown 
that other interaction terms are formally of higher order than the calculated ones. These 
calculations can be systematically improved in the future, if desired. 

In the second part of the paper, we have performed detailed numerical studies of 
the effective theory for the cases of SU{2) and SU(3), focusing on the description of 
its order/disorder phase transition, which is related to the deconfmement transition of 
the 4d thermal theories. We computed the critical couplings and found the transition to 
correspond to 3d Ising for SU(2) and first order for SU(3). Explicit numerical checks with 
different approximations show that next-to-nearest neighbour couplings have negligible 
numerical effects, and furthermore demonstrate good convergence behaviour of the strong 
coupling series. Using the analytic expressions calculated earlier, the critical couplings can 
be converted to values for f3 c (N T ) predicting the deconfmement transition in the original 
4d thermal theories. For lattices with N T = 4 — 16, we found no more than 6% deviations 
from those calculated in simulations of the 4d theories. 

In conclusion, we have given a successful description of the deconfmement transition of 
4d SU(N) Yang-Mills theories in terms of a dimensionally reduced Z (N)-model derived by 
a strong coupling expansions. Future work might push for increased precision by inclusion 
of higher orders, which would entail more interaction terms and higher representation loops. 
Our model might also be of interest for numerical investigations of the deconfmement 
transition for SU(N) with N > 3, cf. [36] and references therein. However, we believe 
at this stage it would be most interesting to extend this approach towards physical QCD 
by including fermions and finite baryon density, e.g. by means of a hopping parameter 
expansion [25, 37]. 
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Corrected integral measure for 577(3) 



The SU(3) effective theories derived by strong coupling expansions in [38] have a path 
integral representation with a measure given in terms of traced Polyakov loops, Li = L(xi), 

Z = (II / dL -) e ~ ScS > s es = ~ Yl + ^AiReLjL*) -^V x . (5.1) 

The potential term in S e g is the Jacobian induced by the Haar measure of the original 
group integration. For our parametrisation L(9,(f)), eq. (3.2) in [38], eq. (3.4) of [38] is 
incorrect. The corrected formula reads 

dL x e v * = f + d6 x [ + d<fe e 2V * , (5.2) 

J —TV J —TV 

with V x given in eq. (3.3); see also eq. (18) in [18]. 

The corresponding corrected effective theories investigated numerically are those with 
(1) one-coupling, (1,2) nearest- and next-to-nearest-neighbour couplings, (l,a) fundamen- 
tal- and adjoint-nearest-neighbours (summations on <, > are on nearest-neighbours, those 
on [, ] are on next-to- nearest neighbours): 

z w = nf d9 *J d ^ n ( 2? - i8 ' l -i 2 + 8Re ( L ') - i l -i 4 ) n ^ + ^iRe^L*) 



<l,3> 



swa n»«p(2Vi,) 

Z{x,2) = [ V[6,0\l[eM2Vx) II (1 + 2AiReL i L*) J](l + 2A 2 ReL fe L;) 

x <i,j> [k,l] 

Z(i,a) = fv[e,4>]Y[exp(2V x ) n (1 + 2A 1 ReL i L*) ]J [1 + \ a {\L m \ 2 - \){\L n \ 2 - 1)] . 

x <i,j> <m,n> 

Clearly, this change in the effective action affects all numerical results obtained in [38] 
for the case of SU(3). It turns out that the correction actually facilitates the simulations 
while modifying the results at the few percent level, such that the conclusions given in [38] 
remain the same. However, for completeness we provide corrected plots obtained from the 
actions above in the following section (see also [39]). All formulae and numerical results 
for the case of ST/ (2) in [38] are correct. 



Corrected numerical results 

In section (3.2) of [38] we discussed the occurrence of a sign problem due to the factor (1 + 
2AiReTjL*), and proposed a solution by approximating the corresponding exponentiated 
logarithm, explog(l + 2AiReLjLp, by a few terms of its series expansion up to order M, 
eq. (3.6) in [38]. Fortunately, with the corrected measure the average sign of this factor 
remains > 0.999 for the couplings of interest, < A < A c ~ 0.188, on all volumes simulated. 
This is demonstrated in figure 17. Hence, the corrected actions can be simulated exactly, 
such that section 3.2 of [38] as well as approximations of the logarithm by polynomials of 
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Figure 17. Average sign of (1 + 2XiReLiL*). Left: a wide range of Ai on a N s 
zoom on the transition region with different system sizes. 
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Fte(L) 

Figure 18. Left: Distribution of L for small and large A*. Middle, Right: Expectation value of 
\L\ and its susceptibility. The vertical line marks the infinite- volume transition. Replaces figure 8 
in [38]. 

order M are now obsolete. In the sequel we list updated plots and tables replacing the 
previous ones. We also report the corrected result for the one-coupling critical value (see 
figure 19): 

\M = 0.187885(30) . (5.3) 

Note that this value differs from the one obtained for the M = 1-truncation of the theory, 
A^ M_1 ^ = 0.13721(5) [40]. The corrected critical lines for the two-coupling theories are 
shown in figure 20. The j3 c {N T ) of the original theory are based on the maps eqs. (2.21-2.23) 
in [38] and eq. (5) in [39], with corrected values given in table 1. 
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